Relativistic effects in two valence electron atoms and ions and search for variation of 

the fine structure constant 
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We perform accurate calculations of the dependence of transition frequencies in two valence elec- 
tron atoms and ions on a variation of the fine structure constant, a — e 2 /he. The relativistic 
Hartree-Fock method is used with many-body perturbation theory and configuration interaction 
methods to calculate transition frequencies. The results are to be used in atomic-clock-type labora- 
tory experiments designed to test whether a varies in time. 
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I. INTRODUCTION 

Theories unifying gravity with other interactions allow 
for the possible variation of physical constants (see, e.g. 
0,0, E|). Recent analysis of quasar absorption spectra 
suggests that the fine structure constant a might vary 
in space-time [3-00. There is an intensive search for 
alternative ways to test whether a is varying. One of 
the very promising methods to study local present-day 
variation of fundamental constants in time involves the 
use of atomic clocks. In particular, optical atomic clock 
transitions are suitable to study the possible variation 
of the fine structure constant. This is because the ratio 
of the frequencies of the optical transitions depend on a 
alone, while the frequencies of the hyperfine transitions 
also depend on the nuclear magnetic moments and the 
electron-proton mass ratio. 

Laboratory measurements involve measuring how the 
difference between two frequencies changes with time. To 
relate a measurement of the change between two frequen- 
cies to a change in a, the relativistic energy shifts are 
needed. The relativistic energy shift describes how a level 
moves as a varies. Two transition frequencies with very 
different relativistic energy shifts are the most desirable 
candidates for precision experiments as they will have the 
largest relative frequency shift between them. 

The best limit on local present day variation of the 
fine structure constant published to date was obtained by 
comparing cesium and rubidium atomic fountain clocks 
Q. Experiments have also been carried out compar- 
ing cesium and magnesium Q and a H-maser compared 
with a Hg II clock There are many proposals for 
the search of variation of a in atomic optical transitions, 
some of which were analyzed previously in 0, 0, ^| . 
In the present work we perform relativistic many-body 
calculations to find the relativistic energy shift for many 
two valence electron atoms and ions. Two valence elec- 
tron atoms and ions were chosen since many new optical 



clocks experiments, some of which are currently under 
construction and some still under consideration, utilize 
these atoms and ions (e. g.. Al II [lg, Ca I UM . Sr I 

EEEzL in n EEIEI, Yb i, H g i HJifT 



II. THEORY 

In the present work we perform calculations for closed 
shell atoms and ions which can also be considered 
as atoms/ions with two valence electrons above closed 
shells. We start our calculations from the relativistic 
Hartree-Fock (RHF) (also known as Dirac-Hartree-Fock) 
method in the V N approximation. This means that RHF 
calculations are done for the ground state of the corre- 
sponding atom/ion with all electrons included in the self- 
consistent field. The use of the V N RHF approximation 
ensures good convergence of the consequent configura- 
tion interaction (CI) calculations for the ground state. 
Good accuracy for excited states is achieved by using a 
large set of single-electron states. Note that there is an 
alternative approach which uses the V N ~ 2 starting ap- 
proximation (with two valence electrons removed from 
the RHF calculations). This approach has some advan- 
tages, it is simpler, and ground and excited states are 
treated equally. However, the convergence with respect 
to the size of the basis is not as good and the final results 
are better in the V N approximation. We use the V N ~ 2 
approximation as a test of the accuracy of calculations of 
the relativistic energy shifts, while presenting all results 
in the V N approximation. 

We use a form of the single-electron wave function that 
explicitly includes a dependence on a: 



V(r) 



_ 1 { f(r) n tt(r/r)ji m 
r V iag(r) n Q{r/r)ji r 



(1) 



This leads to the following form of the RHF equations 
(in atomic units): 
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/n(r) + -ffn(r) - [2 + a 2 (e n - V HF )]g n {r) = 0, 

9n(r) + —g n (r) + (e n - V HF )f n (r) = 0, (2) 
r 
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where k = (— l)' +3 ' +1 / 2 (j + 1/2), n is the principle quan- 
tum number and Vhf is the Hartree-Fock potential. The 
non-relativistic limit corresponds to setting a = 0. 

We then use the combination of the configuration in- 
teraction (CI) method with the many-body perturbation 
theory (MBPT)[2iJ 0|- Interactions between valence 
electrons are treated using the CI method while correla- 
tions between the valence electrons and the core electrons 
are included by means of the MBPT. We can write the 
effective CI Hamiltonian for two valence electrons as: 

H CI = hi + h-i + h u (3) 

here hi (i = 1 or 2) is an effective single-electron Hamil- 
tonian given by 

Ze 2 

hi=caxp + (0- l)mc 2 h V core + E x , (4) 

n 

Vcore is the Hartree-Fock potential created by the core 
electrons, it differs from Vhf in Eq. (J2J by the contribu- 
tion of the valence electrons. Si is the one-electron op- 
erator that describes the correlation interaction between 
a valence electron and the core. The third term in Eq. 
describes the interaction of the valence electrons with 
each other and can be written as 

e 2 

h 12 = — + E 2 (5) 
r\2 

where £ 2 is a two-particle operator that describes the ef- 
fects of screening of the Coulomb interaction between the 
valence electrons by the core electrons. The operators £1 
and E 2 are calculated using the second order of MBPT. 

We use the same set of single-electron basis states to 
construct two-electron wave functions for the CI calcula- 
tions and to calculate £. The set is based on the B-spline 
technique developed by Johnson et al [2^, |2|| |53 • We 
use 40 B-splines in a cavity of radius R = AQcib (as is 
Bohr radius). The single-electron basis functions are lin- 
ear combinations of 40 B-splines and are also eigenstates 
of the Hartree-Fock Hamiltonian (in the V N potential) . 
Therefore, we have 40 basis functions in each partial wave 
including the B-spline approximations to the atomic core 
states. We use a different number of basis states for the 
CI wave functions and for the calculations of £. Satura- 
tion comes much faster for the CI calculations. In these 
calculations we use 14 states above the core in each par- 
tial wave up to l max — 3. Inclusion of states of higher 
principal quantum number or angular momentum does 
not change the result. To calculate E we use 30 out of 
40 states in each partial wave up to l ma x = 4. 

The results for the energies are presented in Table [I] 
We present the energies of the nsnp configuration of two 
electron atoms/ions with respect to their ground state 
1 S'o ns 2 . The states considered for atomic clock experi- 
ments are 3 Po an d 3 Pi- However, we present the result 
for other states as well for completeness, these also make 
it easier to analyze the accuracy of the calculations. Also, 



TABLE I: Energies of the nsnp configuration of two electron 
atoms calculated using H CI , H CI + Ei and H CI + Ei + E 2 ; 
comparison with experiment (cm -1 ) 



Atom/ State Experiment Theory 



ion 




[30] 


jjCI 


H CI + Ei 


H CI + E 


A1II 


6 Po 


37393 


36403 


36987 


37328 




3 Pi 


37454 


36466 


37053 


37393 




3 P 2 


37578 


36592 


37185 


37524 




'Pi 


59852 


59794 


60647 


60090 


Cal 


3 Po 


15158 


13701 


14823 


15011 




3 P1 


15210 


13750 


14881 


15066 




3 P 2 


15316 


13851 


14997 


15179 




'Pi 


23652 


23212 


24968 


24378 


Sri 


3 Po 


14318 


12489 


13897 


14169 




3 P 


14504 


12661 


14107 


14367 




3 P 2 


14899 


13021 


14545 


14786 




'Pi 


21698 


20833 


23012 


22305 


Inll 


3 Po 


42276 


37825 


39238 


42304 




3 P 


43349 


38867 


40394 


43383 




3 P 2 


45827 


41168 


42974 


45904 




'Pi 


63034 


62181 


64930 


62325 


Ybl 


3 Po 


17288 


14377 


16352 


16950 




3 P 


17992 


15039 


17189 


17705 




3 P 2 


19710 


16550 


19137 


19553 




'Pi 


25068 


24231 


27413 


26654 


Hgl 


3 Po 


37645 


31864 


32692 


37420 




3 P 


39412 


33751 


34778 


39299 




3 P 2 


44043 


38155 


39781 


44158 




'Pi 


54069 


50247 


52994 


56219 


Till 


3 Po 


49451 


43831 


43911 


49865 




3 P 


52393 


47091 


47350 


52687 




3 P 2 


61725 


55988 


56891 


62263 




'Pi 


75660 


74291 


76049 


74717 



transitions associated with some of these states are ob- 
served in quasar absorption spectra (e.g., the So — 1 Pi 
transition in Ca). 

To demonstrate the importance of the core- valence cor- 
relations we include results of pure CI calculations (with 
no E) as well as the results in which only £1 is included 
but E2 is not. One can see that the accuracy of pure CI 
calculations is about 10% while inclusion of core-valence 
correlations improves it significantly to the level of about 
1%. The deviation from experiment of the final theoreti- 
cal energies for the triplet states of all atoms except Yb is 
not more than 1%. For Yb it is 2%. The accuracy of the 
singlet states is about 1% for the ions, 3-4% for Cal, Sri 
and Hgl and 6% for Ybl. The accuracy of the fine struc- 
ture intervals ranges from 2 to 7%. The accuracy of cal- 
culations for Yb is not as good as for other atoms because 
the two electron approximation is a poor approximation 
for this atom. Electrons from the 4/ subshell, which are 
kept frozen in present calculations, are relatively easy to 
excite and corresponding configurations give substantial 
contribution to the energy. Note that we do include these 
excitations perturbatively, into the E operator. However, 
due to their large contribution, second-order treatment 
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TABLE II: Calculated q coefficients, for transitions from the 
ground state, using H CI , H CI + Ei and H CI + Ei + E 2 



Atom/ion 


State 




H + Ei 


H°' + Ei, 2 


Other 


A1II 


3 Po 


138 


142 


146 






3 Pi 


200 


207 


211 






3 P 2 


325 


340 


343 






1 Pl 


266 


276 


278 




Cal 


3 Po 


108 


115 


125 






3 Pi 


158 


173 


180 


230 [10] 




3 P 2 


260 


291 


294 




*Pi 


228 


238 


250 


300 [10] 


Sri 


3 7 1 

3 Po 


384 


396 


443 




3 Pi 


560 


609 


642 


667 [31] 




7 — . 

3 P 2 


939 


1072 


1084 




X Pl 


834 


865 


924 


1058 [31] 


Inll 


3 Po 


3230 


2932 


3787 


4414 [12] 




^ 7-1 

3 Pi 


4325 


4125 


4860 


5323 [12] 




3 P 2 


6976 


7066 


7767 


7801 [12] 




1 7 1 

X Pl 


6147 


6103 


6467 


Ybl 


'i 7 1 

3 Po 


2339 


2299 


2714 






3 Pi 


3076 


3238 


3527 






3 P 2 


4935 


5707 


5883 






X Pl 


4176 


4674 


4951 




Hgl 


3 Po 


13231 


9513 


15299 






3 Pl 


15922 


12167 


17584 






3 P 2 


22994 


19515 


24908 






J Pl 


20536 


16622 


22789 




Till 


3 Po 


14535 


11101 


16267 


19745 [12] 




3 P1 


18476 


14955 


18845 


23213 [12] 




3 P 2 


32287 


28903 


33268 


31645 [12] 




X P1 


28681 


25160 


29418 



TABLE III: Experimental energies and calculated q- 
coefficients (cm -1 ) for transitions from the ground state ns 2 
to the nsnp configurations of two-electron atoms/ions 



Atom/Ion 


Z 


State 


Energy [30] 


q 


A1II 


13 


3s3p 


3 Po 


37393.03 


146 






3s3p 


3 Pi 


37453.91 


211 






3s3p 


3 P 2 


37577.79 


343 






3s3p 


J Pi 


59852.02 


278 


Cal 


20 


4s4p 


3 Po 


15157.90 


125 






4s4p 


3 P1 


15210.06 


180 






4s4p 


3 P 2 


15315.94 


294 






4s4p 


X P1 


23652.30 


250 


Sri 


38 


5s5p 


3 Po 


14317.52 


443 






5s5p 


3 P1 


14504.35 


642 






5s5p 


3 P 2 


14898.56 


1084 






5s5p 


'Pi 


21698.48 


924 


Inll 


49 


5s5p 


3 Po 


42275 


3787 






5s5p 


3 P1 


43349 


4860 






5s5p 


3 P 2 


45827 


7767 






5s5p 


'Pi 


63033.81 


6467 


Ybl 


70 


6s6p 


3 Po 


17288.44 


2714 






6s6p 


3 P1 


17992.01 


3527 






6s6p 


3 P 2 


19710.39 


5883 






6s6p 


'Pi 


25068.22 


4951 


Hgl 


80 


6s6p 


3 Po 


37645.08 


15299 






6s6p 


3 P1 


39412.30 


17584 






6s6p 


3 P 2 


44042.98 


24908 






6s6p 


'Pi 


54068.78 


22789 


Till 


81 


6s6p 


3 Po 


49451 


16267 






6s6p 


3 P1 


53393 


18845 






6s6p 


3 P 2 


61725 


33268 






6s6p 


'Pi 


75600 


29418 



of the excitations from the 4/ subshell is not very accu- 
rate. On the other hand, the CI+MBPT results for Yb 
are still much better than pure CI values. 

Note also that the CI+MBPT results presented in Ta- 
ble [I] are in g ood agreement with similar calculations in 
Refs. [HHa]. 

III. RESULTS AND DISCUSSION 

In the vicinity of the ao, the present day value of a, 
the frequency of a transition, u, can be written as: 

lu = ujq + qx, (6) 

where x = (-^-) 2 ~ 1, u>q is the present day experimental 
value of the frequency and the q coefficient is the rela- 
tivistic energy shift that determines the frequency depen- 
dence on a. It is clear from the above expression that q 
coefficients can be described by 

duj . 

q = -7- \x=o- 

ax 

Thus, in order to calculate q coefficients the atomic en- 
ergy levels of the atoms and ions of interest at different 



values of x need to be calculated. The relativistic energy 
shift q is then calculated using the formulae 

q = u(Ax) -u(-Ax) ^ 

and 

_ l6(w(Ax) - uj(-Ax)) - 2{oj(2Ax) - oj(-2Ax)) 

q ~ 24Ax ■ 

(8) 

The second formula is needed to check for non-linear con- 
tributions to dui/dx. We use Ax = 0.1 and Ax — 0.125. 
The results are presented in Table HP 

As for the energies, we use three different approxima- 
tions to calculate relativistic energy shifts: (1) pure CI 
approximation for two valence electrons, (2) CI with Si 
and (3) CI+MBPT approximation with both Ei and £2 
included. Inclusion of core- valence correlations lead to 
increased values of the g-coefficients. This is because 
the correlation interaction of a valence electron with the 
core introduces an additional attraction which increase 
the density of the valence electron in the vicinity of the 
nucleus and thus emphasize the importance of the rela- 
tivistic effects. 
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Note that Si and £2 are of the same order and need 
to be included simultaneously to obtain reliable results. 
Ei is much easier to calculate and inclusion of £1 alone 
often leads to significant improvements of the results for 
the energies (see Table HJ. However, the results for the 
q-coefficients show that neglecting £2 may lead to sig- 
nificant loss in accuracy. Indeed, the results for g's with 
£1 alone are often smaller than those obtained in pure 
CI and CI+MBPT approximations and differ from final 
values by up to 50%. Since neglecting £2 cannot be jus- 
tified, we present results without £2 for illustration pur- 
pose only. 

The accuracy of the calculation of the ^-coefficients can 
be estimated by comparing the CI and CI+MBPT results 
calculated in the V N and V N ~ 2 approximations and also 
by comparing the final results for the energies (including 
fine structure intervals) with experimental values. As 
one can see from Table [n] inclusion of the core- valence 
correlations can change the values of the q-coefficients 
by more than 15%. However, the accuracy of the ener- 
gies improves significantly when core- valence correlations 
are included. It is natural to expect that the final accu- 



racy for the ^-coefficients is also higher when core- valence 
correlations are included. Comparison with our previ- 
ous results also shows some deviation on approximately 
the same level (the largest relative discrepancy is for Ca 
where relativistic effects are small and high accuracy is 
not needed). Most of this discrepancy can be attributed 
to the inaccuracy of our old, less complete calculations. 
Comparison between the energies calculated in the V N 
and V N ~ 2 approximations and the experimental values 
suggest that 10% is a reasonable estimate of the accuracy 
of the present calculations of the relativistic energy shifts 
for Al II, Ca I and Sr I, 15% for In II, 25% for Yb I and 
20% for Hg I and Tl II. 

In Table ITTT1 we present final values of the relativistic 
energy shifts together with the experimental energies. 
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